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Abstract 

We discuss a generalization of the minimal supersymmetric extension of the Stan- 
dard Model in the form of three additional singlet superfields, which would explain 
the essential features of the CDF multi-muon events presented recently: a large pro- 
duction cross section of ~ 100 pb originates from the production of a CP-odd scalar 
A with a mass in the 70 - 80 GeV range and a large value of tan f3 ~ 40. The CP-odd 
scalar A decays dominantly into CP-odd and CP-even scalars ai and hi, which gener- 
ate decay cascades hi 2/i2 — > 4o2 8t, and oi — > hia2 with hi decaying as above. 
The decay 02 — > r"*"T~ is slow, leading to a lifetime of 0(20) ps. The phenomenology 
of the model differs from similar scenarios presented before in that one of the two 
cascades leads to 10 instead of 8 r-leptons, and additional production processes like 
associate A production with bb pairs are relevant. 
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1 Introduction 



In a recent CDF publication [1] a sample of unusual multi-muon events was studied, which 
are characterized by the following properties: i) a large rate of additional muons inside 
a cone of 36.8° (cos 6' > 0.8) around the direction of the trigger muon; ii) large impact 
parameters (displaced vertices), and iii) an unusually large cross section of the order of 
75 pbQ Since neither Standard Model processes nor known detector effects can explain the 
nature of these events at present, they were refered to as ghost events. 

In [2], some members of the CDF collaboration investigated to which extent the prop- 
erties of the ghost events can be understood in the context of a phenomenological scenario 
based on cascade decays of new particles. First, the sign-coded multiplicity distribution of 
additional muons inside the cos^ > 0.8 cone coincides with the assumption that originally 
4r"'" + 4r~ leptons were produced. As a hypothetical origin for the 8 r-leptons, the authors 
of [2] considered the pair production of new particles hi via pp ^ H hihi, without 
specifying the nature of H (which could be a known or a new gauge boson, or another new 
particle). Subsequently, each of the hi particles is assumed to produce the decay cascade 
hi 2/i2 — *• 4/i3 —> 4(r"'" + r~), generating two multi-muon cones per event. /i2 and /is 
denote additional new states, with h^ decaying as /13 — *■ r"*" -|- r~. 

In order to explain the high multiplicity of additional muons inside the cos 6* > 0.8 
cones, the particle hi must be relatively light. The best fit [2] to the invariant mass 
distributions inside the cos^ > 0.8 cones in [1] is obtained for hi, h2 and h^ masses near 
the lower limits for which the cascade is kinematically allowed: m/13 ~ 3.6 GeV > 2mT-, 
mfi2 ~ 7.3 GeV > 2m/i3 ''^hi ~ 15 GeV > 2mh^. In order to explain the large impact 
parameters, at least one of the hi, h2 or /13 particles must have a long lifetime. The best 
fit [2] to the impact parameter distributions in [1] corresponds to the assumption that h^ 
has a long lifetime of ~ 20 ps. The origin of the large cross section (for the production of 
H in this scenario) was left unexplained in [2]. 

Some comments on the multi-muon study by the CDF collaboration [1] and the phe- 
nomenological interpretation in [2] were published in [3], wherein proposals for additional 
studies/plots were made, and some difficulties with the phenomenological scenario in [2] 
were pointed out, on which we will comment later. Furthermore, different phenomenolog- 
ical scenarios (various types of micro-cascades hi ff'hi^i) were proposed in [3], which 
could provide a better fit to the data. 

Theoretical models, in which multi-lepton events at hadron colliders can be expected, 
have already been considered before [4-8]. One class of such models contains a rich nearly 
hidden sector (a "hidden valley"), which interacts only weakly with the Standard Model 
(SM) particles, as via a heavy Z' gauge boson, or via a small kinetic mixing between the SM 
and the hidden sector gauge fields [4-8]. It seems possible to explain the recently observed 
high-energy components of cosmic rays as remnants of dark matter annihilation in such 
models [9]. (See [10] for a discussion of string vacua including D-branes, where a light 
hyperweak gauge boson can connect the SM with a hidden sector, and ghost-like events 
can occur.) 

^The precise value of the necessary total cross section depends on the fraction of the events which survive 
the cuts applied in [1], which requires a model specific simulation. 
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Multi-lepton events at hadron colliders are also possible within the Next-to-Minimal 
Supersymmetric Standard Model (NMSSM) [11-17], where the Higgs sector of the MSSM 
is extended by a singlet superfield S. The extended Higgs sector can contain a light CP- 
odd state Oi with a mass below the bb threshold of ~ 10.5 GeV, such that ai decays 
dominantly as Oi t^t' . If, in addition, the SM-like Higgs scalar hsM decays dominantly 
as hsM ^ fli + Oi, events with 4 r-leptons could be the only signal of hsM^ rendering 
its discovery quite difficult (see [18-20] and refs. therein). Alternatively, two leptons can 
originate from a bino decay into a singlino-like Lightest Supersymmetric Particle (LSP) [21, 
22]. However, in both cases the large number of muons observed by CDF is not achieved. 

In the present paper we consider the extension of the Higgs sector of the MSSM by 
several (three) singlets in order to discuss whether the properties of the CDF ghost events 
could be understood in such a setup. An extended Higgs sector was already implicitely 
suggested by the authors of [2] as a source of the states /ii, /i2 and /13, and the construction 
of corresponding models leading to the desired masses in the range 3.5 - 15 GeV does 
not seem very difficult at first sight. However, in practice various problems appear: First, 
the unusually large production cross section must be explained. Second, the desired decay 
channels must be dominant, and at least one large lifetime must occur. Third, and most 
importantly, present constraints from colliders (notably LEP), B physics etc. on such an 
extended Higgs sector must be satisfied. In the worst case it may be impossible to satisfy 
all these conditions simultaneously, at least once the masses and couplings of the extended 
Higgs sector are constrained by supersymmetry. 

Hence, it is important to look for a "go-theorem" : a concrete model, which has all these 
desired properties. Here we present such a model, which has the following structure: The 
starting point is the NMSSM involving a singlet superfield 5, whose vacuum expectation 
value (vev) solves the //-problem of the MSSM. We consider a region in the parameter space 
of the NMSSM, where the lightest CP-odd Higgs state in the Hu-Hd-S sector (denoted by 
A subsequently) has a mass in the 70 - 80 GeV range and has both large singlet and large 
doublet components of ~ 85% and ~ 50% respectively. Assuming a large value of tan j3 ~ 
40, the production cross section of A via gluon-gluon fusion at the Tevatron can be ~ 100 pb 
[23]. In the corresponding region in parameter space, the CP-even Higgs states in the Hu- 
Hd-S sector as well as the second CP-odd state have masses above 114 GeV, in which case 
all LEP constraints on various Higgs production processes [24] are satisfied. (The precise 
Higgs masses as well as various B-physics observables depend on the soft supersymmetry 
breaking gaugino, squark and slepton masses and couplings, for which ranges of desired 
values can be found without particular effort.) 

To the Higgs sector of the NMSSM we add two more singlets 5*1 and 5*2, which contain 
two more CP-even states hi^2 and two more CP-odd states ai.2. (From here onwards, the 
indices 1,2 denote the additional singlets of the model rather than the states introduced 
in [2].) These supplementary fields allow for many additional Yukawa couplings and soft 
terms, but for simplicity we assume that most of the possible Yukawa couplings vanish or 
are negligibly small (considering small Yukawa couplings, as those appearing in the SM, 
as natural). Then, the following situation can be achieved without fine tuning, assuming 
corresponding values of the additional soft terms: The mass matrices of the hi 2 and 01^2 
states are nearly diagonal, with eigenvalues in the 3.5 - 20 GeV range. A Yukawa coupling 
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between 5* and 5*1 remains relatively large, and as a consequence the CP-odd state A of the 
Hu-Hd-S sector decays dominantly into hi + ai (as compared to A ^ bb). Subsequently, 
due to a small S1-S2 mixing, hi decays as hi 2/i2 — > 4a2 — >■ 8r, and ai as ai — > hia2 with 
/ii-decays as above. 

The decay 02 t^t" is possible due to a tiny S2 — mixing of C(10~^), implying a 
a2 lifetime of C(20) ps. Assuming m^^ ~ 20 GeV and nih-^ ~ 15 GeV, the first A hiai 
decay will generate two separate cones containing 8 or 10 r-leptons (which are thus not 
completely symmetric) as well as displaced vertices. Thus, the essential properties of the 
phenomenological scenario of [2] are reproduced, and the required cross section is obtained. 

On the other hand the concrete model makes it clear that a phenomenological analysis 
based on a single process can be incomplete or even misleading: In the present scenario, 
the CP-odd state A will also be produced in association with bb pairs with a cross section 
of ~ 30% of the production via gluon-gluon fusion [23]. In ~ 25% of these cases, a 6 or 
b decay will generate at least one additional muon, which can end up in one of the cones 
defined by the trigger muons. These additional muons will have an impact on observables 
like \Pt\ and invariant masses, whose precise effect can only be studied with the help of 
simulations, which are beyond the scope of the present paper. (Likewise, the production and 
the decays of the heavier states of the H^-Hd-S sector can contribute to the observables.) 
As already underlined in [3] , an analysis of the data involving more stringent cuts than the 
ones presented in [1] (as muon charge selection rules) should make it easier to constrain - 
or to verify - concrete models as the one discussed here. 

We are aware of the fact that the properties of the CDF multi-muon event samples 
still need to be confirmed, notably by the DO collaboration. Nevertheless we found it 
useful to develop a concrete model, which indicates which additional complications can be 
expected and which, in any case, extends the perimeter of possible signals that may be 
expected at hadron colliders. Of course, if the properties of the multi-muon event samples 
are confirmed, it is also interesting to know that relatively simple singlet extensions of the 
MSSM can generate such signals. In the next section we will present the Lagrangian and 
discuss the parameters of the model. Its phenomenology and conclusions will be presented 
in section 3. 

2 A toy model 

As mentioned in the introduction, we consider a supersymmetric extension of the SM in- 
volving an extended Higgs sector, which consists of the MSSM doublets Hu, Hd and three 
singlet superfields S", Si and 82- A priori, a large number of terms could appear in the su- 
perpotential W , even after the restriction to scale invariant Yukawa couplings as motivated 
by a solution of the /i-problem of the MSSM. On the other hand, small or vanishing Yukawa 
couplings are "technically natural" (stable under quantum corrections) in supersymmetry, 
and we use this freedom to omit most of the allowed couplings. 

The relevant Higgs dependent terms in the superpotential W are assumed to be given 

by 

W = XSH^Hd + + XiSS! + '^S! + X2SiSi + '^Sl . (2.1) 
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The first two terms correspond to the ones of the NMSSM [11-17], where the vev s = (S) 
generates an effective /i-parameter fics = As. The terms proportional to ki and K2 serve to 
stabilize the scalar potential for the vevs Si = (Si) and S2 = {S2), and Ai and A2 induce 
couplings between S and 5*1, and between 5*1 and S2, respectively. 

The soft terms in the Higgs sector are the scalar masses squared and trilinear couplings 

^Hu^ "^L' "^55 "^5i> ^S2^ ^A, A^, Ax,, A^-,, Ax^ and v4«2 . (2.2) 

It is straightforward to work out the tree level Higgs mass matrices and couplings from 
the superpotential and the soft terms, which leads to quite lengthy and not very transparent 
expressions. Instead of presenting them, we first consider the "decoupling limit" A2 — >■ 0. As 
it becomes evident from the superpotential, all components of the superfield 5*2 decouple 
from S and Si in this limit. Furthermore one finds that for a wide range of 

parameters (for A\^ not too large and positive m\ , see eq. (12.31) below), the vev si vanishes, 
since terms linear in si in the scalar potential are proportional to A2. Then the mass matrices 
in the Higgs sector are block- diagonal: 

In the NMSSM sector Hu-Hd-S one re-obtains the well known 3x3 (2x2) mass matrices 
for the CP-even (CP-odd) states [11-17]. In addition, the mass matrices for the CP-even 
hi^2 and CP-odd 01,2 states are diagonal. In the case of the vev S2, we assume that \Ai^^\ is 
sufficiently large such that the vev $2 is nonvanishing (which avoids degenerate /i2, a2 states 
and a massless neutralino il)2)- Then it is convenient to express m^^ in terms of S2, ^2 and 
through the minimization equation of the scalar potential, after which the masses of 
the physical states hi^2 and 01^2 can be written as 

= ml^+2XiAx,s + 2KXis^ -2XXiVuVd + 4:Xls^ , 
m1_^ = — 2XiAx,s — 2kXis^ + 2XXiVuVd + 4A^s^ , 

"^h2 = '«2S2(^K2 + 4K2S2) , 

"^aa = -3fi:2v4^2S2 (2.3) 

where Vu-, Vd denote the vevs of the neutral components of H^, Hd- 

Evidently there exist sufficient free parameters Ai, Ax,., S2, ^2 and in the 5'i-S'2 
sector in order to generate a spectrum like 

ma, ~ 20 GeV , nih, ~ 16 GeV , m/,^ ~ 8 GeV , m^^ ~ 4 GeV , (2.4) 

which render the cascade decays described in the introduction kinematically possible, with 
relatively light initial states hi and ai. (The masses of the additional neutralinos are given 
by = 2Ais and m^^ = 2^282- "^i is too heavy to be produced in A decays, and the 
BR(yl — i> ip2'ip2) vanishes in the decoupling limit.) 

Of course, the desired cascade decays of hi and ai require the presence of couplings 
9hih2h2^ 9aihia2 ga2T+T-, which are absent in the decoupling limit A2 — * 0. (The coupling 
9h2a2a2 is of the order i^2Ak2 and not suppressed in the decoupling limit.) One can check 
that a small value of A2 will generate couplings of the order (modulo Yukawa couplings and 
a dimensionful parameter like an A-term or a vev) (7/11/12/12 ~ 9aihia2 ~ A2, QaiT+r- ~ A|. 
One finds that for Ag ~ 10~^, the 02 lifetime will be of the order 20 ps as desired, but which 
has no noticeable effect on the eigenvalues of the mass matrices above. 
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In the NMSSM sector, the tree level mass matrices receive considerable radiative correc- 
tions depending on the squark, slepton and gaugino masses and trilinear couplings. These 
are included in the code NMHDECAY/NMSSMTools [25-27], which computes the NMSSM 
Higgs masses and couplings as functions of the parameters in the Lagrangian. As indepen- 
dent parameters in the Higgs sector of the NMSSM, one can choose [25-27] 



Large cross sections of Higgs particles at hadron colliders occur at large values of tan /3, 
for which the coupling of Hd to down quarks is proportional to tan (3. Then, the h quark 
loop induced gluon-gluon fusion process has a cross section amplified by ~ tan^/3 with 
respect to the corresponding cross section for the production of a SM Higgs scalar. On the 
other hand, the cross section via gluon-gluon fusion decreases strongly with increasing Higgs 
masses, but low CP-even Higgs masses are strongly constrained by LEP [24]. Therefore 
we concentrate on a region in the NMSSM parameter space (12. 5p at large tan/5 where a 
CP-odd Higgs scalar A has a mass below 100 GeV, but large enough to render the 
decay A — > hiai (with A on-shell) possible. Also, the decay A —>■ hiai should have a larger 
branching ratio than the decay A —>■ bb, which requires a considerable singlet 5* component 
of A without a too large suppression of the coupling of A to 6 quarks (see the next section). 

Finally, LEP constraints on all CP-even Higgs scalars as well as constraints from 
B physics (the branching ratios BR(i? — s> Xg'y), BR(i?"^ — > T~^i>r), BIi{Bs f^~^f^~) and the 
mass differences AMg, q = d,s), which are particularly relevant at large tan (3, should be 
satisfied. All these constraints are checked in the code NMHDECAY/NMSSMTools [25- 
27], which we used in the search for acceptable regions in the parameter space (12.51) of 
the model. Clearly, Higgs masses as well as B physics observables depend also on the soft 
supersymmetry breaking gaugino, squark and slepton masses and couplings, which have to 
be specified. 

In fact, for large tan/5 ~ 40 and non-negligible NMSSM Yukawa couplings A and k a 
phenomenologically acceptable region in the parameter space exists, in which A has a mass 
in the 70 - 80 GeV range and Hd- and ^-components of ~ 50% and ~ 85%, respectively. 
An example is given by the following point in parameter space, where 

A = .28, K = .33, A^ = -36.5, tan/5 = 40, /ieg = 240 GeV, Ma = 420 GeV . (2.6) 

The gaugino masses are Mi = 150 GeV, = 300 GeV, M3 = 1 TeV, the left-handed 
and right-handed up- type squark masses are 1.5 TeV, the right-handed down-type squark 
masses are given by 1 TeV, the slepton masses by 500 GeV, Atop = Abottom = 1.8 TeV and 
Ar = 300 GeV. 

For these parameters, the lightest CP-odd Higgs mass niA and its decomposition A = 
Na,Au^u + NA^Ad^d + Na,s^s (where Au, Ad and Ag are the neutral CP-odd components 
of Hu, Hd and S) are given by 



The masses of the three CP-even Higgs scalars are 114.5 GeV, 270 GeV and 561 GeV, 
and the masses of the second CP-odd and charged Higgs scalars are 300 GeV and 264 GeV, 
respectively. Further properties of this point in parameter space, which are relevant for the 
CDF ghost events, will be discussed in the next section. 



A, K, A,^, tan/3, fics = As, = 



2fics{Ax + KS) 

sin 2/? 



(2.5) 



ruA = 70 GeV, Na,a^ = 0.01, Na,a, = 0.56, A^^,^ = 0.83 . 



(2.7) 
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3 Phenomenology of the toy model 



In order to estimate the production cross section of A via gluon-gluon fusion at the Tevatron, 
one has to determine its reduced couphng to down-type quarks (normahzed to the SM 
Higgs couphng), 

Xd = tanPx Na,a, (=22.2), (3.1) 

where the value in parenthesis is the one for the point given in (12.61) . (12.71) . Then, the 
corresponding production cross sections for the Tevatron in [23] can be rescaled appropri- 
ately, and extrapolated to = 70 - 80 GeV. For = 70 GeV and ~ 22 one obtains 
a{pp A + X) ~ 100 pb, even somewhat larger than required. 

Subsequently, we have to estimate the A decay branching fractions. In the absence of 
the 5*1^2 sector, A would decay to ~ 90% into bb with a partial width 



X^mtmAJl (3.2) 



which gets enhanced by ~ 20% by QCD corrections. In the presence of a coupling gAhic 
the partial width for the decay A — ^ hiai is 



r.,a, = -pi^W 1-^-^ -4^V^. (3.3) 
lonmA V V ^aJ ''^a 

Numerically, one obtains for the ratio 



T/iiai / 36 QAhiai 



bb 



XarriA 
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In order to obtain a branching fraction for A hiai larger than ~ 80%, such that the 
production cross section a{pp — > A — > hiai) is larger than ~ 80 pb, we should have R > 
For the values of X^ and above, this can be obtained for QAhiai ^ 86 GeV. In the 
present model, gAhiai is given by 

gAh^ai = (^^1 + ^'^■s) (3.5) 

with Na,s as in (12. 7p . Note that QAhiai is not suppressed in the decoupling limit A2 — > 0. 
It is easy to find values for Ai ~ 0.52 and A\^ ~ —283 GeV such that QAhiai is sufficiently 
large, and the masses mhi and obtained from (12. 3p have the desired values. 

After the dominant decay A hiai, hi and ai cannot decay at tree level in the 
decoupling limit {A2, Si} 0, where hi and ai do not mix with the Hy-H^-S sector. Small 
values of A2 (and appropriate natural values for A\^, k,2 and ^^2) are sufficient in order to 
generate the dominant decays hi — ^ /i2^2; ^2 ~^ fl2'^2 and ai — hia2, which we assume to be 
kinematically allowed, and which produce the cascades described in the introduction. The 



^ We recall the footnote on page 1 . according to which the necessary total cross section can be somewhat 
larger or smaller. 
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decay of 02 into quarks and leptons is made possible only through its small mixing ~ A| with 
the Hu-Hd-S sector. For large tan /3, 02 mixes dominantly with i?^, from which it inherits 
the couplings proportional to the down-type fermion masses leading to the dominant decay 
a2 — »■ T^T~ (for 2mb > rriaj ^ 2mr). Herewith we have reproduced the essential features 
of the scenario proposed in [2], albeit with one of the two cascades (the one originating 
from ai) leading to 10 rather than 8 r-leptons in the final state, which will evidently imply 
some modifications of the plots presented in [2]. 

Next we comment on an issue raised in [3], where it has been noted that a Higgs- 
like coupling to r-leptons implies a coupling to muons with a ratio m^/m^, and hence a 
ratio of branching ratios BR(a2 — ^ /i"'~/i~)/BR(a2 r'^r") > m^/m^ ~ 0.0035. For 
near 2 the ratio of branching ratios increases due to the kinematic suppression of the 
BR(a2 — * r+r"). In plots of the invariant mass of muons of opposite charges inside a cone, 
this could (should) generate a visible peak at ma^- This reasoning remains valid in our 
case, but we note that the kinematic suppression of the BR(a2 — > r+r~) near the threshold 
TTiaj — > 2 rrir is less important for CP-odd scalars (like 02) as compared to CP-even scalars; 
accordingly, the ratio of branching ratios above will increase less dramatically near the 
threshold — > 2m^, which makes it somewhat more difficult to rule out the scenario 
through a non-observation of a peak in the fi~^fi~ invariant mass distribution. 

Finally we turn to the invariant mass distribution M of all muons - or of all tracks - 
for events in which both cones contain at least two muons (Fig. 35 in [1]). According to 
the simulations of the process pp H ^ hihi performed in [2] (in the notation of [2]), a 
resonance-like structure should be visible with a peak position depending on rriH (see Fig. 6 
in [2]); however, the data do not show such a structure: the process simulated in [2] could 
not describe simultaneously the steep rise of the invariant mass distribution for small M, 
and the tail of the invariant mass distribution at large M, for any value of ruH- 

In the present scenario we have to replace H hj A with a mass in the 70 - 80 GeV range, 
which seems to describe only the invariant mass distribution for small M. Furthermore, 
one of the cones would contain 10 r-leptons; however, this is possibly not yet enough in 
order to explain the tail of the invariant mass distribution at large M. 

On the other hand, as already mentioned in the introduction, there exist additional 
production processes which have necessarily to be taken into account: the cross section for 
associate bb + A production can be estimated to be ~ 30% of the A production via gluon- 
gluon fusion (hence ~ 30 pb) for A masses in the range considered here [23]. In ~ 25% of 
these cases, a 6 or 6 decay will generate at least one additional muon which can contribute 
to the tail of the invariant mass distribution at large M. 

Furthermore, one of the heavier CP-even scalars (the one with a mass of ~ 270 GeV for 
the point above) as well as the heavier CP-odd scalar Aheavy (with a mass ~ 300 GeV here) 
have couplings to b quarks enhanced by factors of tan/3 = 40 and tan/? x Na^^^^^^a^ ~ 33, 
respectively. At least in regions in parameter space where their masses are still lower, these 
states - which will generate similar cascades leading to 8 - 10 r-leptons - can also contribute 
to the tail of the invariant mass distribution at large M. Of course, further simulations are 
necessary in order to check these conjectures, but at first sight explanations of the invariant 
mass distributions for small and for large M seem possible. Clearly, these processes will 
also contribute to other observables like YIpt- 
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At last, a comment on the dark matter relic density in this model is appropriate. The 
LSP is the neutralino ip2 with a mass m^^ = 2K2S2 ~ > 2ma^. Its annihilation is 
dominated by the processes ip2 + ^ h2 ^ a2 + 02 and ip2 + 4'2 — >■ /12 + 02- 
Subsequently the scalars 02 (/i2) will decay into two (four) r-leptons as at the end of 
the cascades relevant for the multi-muon events. The annihilation processes depend on 
the 4'2 4'2h2/a2 Yukawa coupling K2, and on the trilinear coupling gh2a2a2 of the order 
of K2v4k2. Whereas A^^ is determined by the desired values of nih^ and fl2.3l) to be 
~ —(1 — 1.5) GeV (and ^2*2 ~ 4 GeV), the value of K2 is unconstrained so far. One 
can expect that, for a value of K2 in the range O{10~^) - O{10~^), the WMAP value 
0.094 < n^^h^ < 0.136 [28] for the dark matter rehc density can be achieved. 

To conclude, apart from the fact that the CDF multi-muon events [1] need to be con- 
firmed notably by the DO collaboration, additional studies of their properties would be 
desirable, as the ones pointed out in [3]: spatial correlations among displaced vertices, and 
invariant mass distributions of dimuon pairs depending on their relative charges. In any 
case, plots have to be compared with simulations of models. 

We have presented a relatively simple model in the form of a multi-singlet extension 
of the MSSM, whose particle content and parameters have been chosen such that the es- 
sential features of the CDF multi-muon events can be reproduced, without contradicting 
constraints from other experiments. Already in this scenario, the phenomenology would be 
more complicated than the one discussed in [2]. Clearly, the particular values of the param- 
eters of the model chosen in section 2 - and the corresponding masses and couplings - have 
been presented for illustrative purposes only, and eventually the complete phenomenologi- 
cally acceptable region in parameter space could be studied. Together with further models, 
which will certainly be proposed soon, this will allow for comparisons or "best fits" to the 
data. 
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